Manifold self-localization in a deformable medium 



> 
o 



o 



I 

O 

o 



o 



On 
ON 



o 
o 



X 
S3 



Eugene B. Kolomeisky^ and Joseph P. Straley^ 
^Department of Physics, University of Virginia, Charlottesville, Virginia 22901 
^Department of Physics and Astronomy, University of Kentucky, Lexington, Kentucky 40506-0055 

Directed manifolds (domain walls, interfaces, vortex lines) in a deformable medium can exist in 
a correlated state in which the manifold is self-localized by its own strain field. Depending on the 
temperature, manifold/medium dimensionalities, and the strength of the couphng with the medium, 
the degree of localization of the ground state can vary both continuously and discontinuously; there 
can be phase transitions from self- localized to the free-manifold state. 



The concept of a manifold is a central paradigm in con- 
densed matter physics, covering the cases of domain walls 
in ferroelectric and ferromagnetic materials, flux lines in 
type-II superconductors, interfaces, and many more [Q. 
Manifolds exist in symmetry-broken states of various sys- 
tems, and occur as translationally-invariant topological 
solutions of the corresponding field theories. 

In this work we show that the deformability of a real 
system can bring the manifold into a new correlated self- 
localized state (SL) which forms through the following 
scenario: a deformable continuum responds to the pres- 
ence of the manifold by a strain field decaying away from 
the manifold; this field interacts with the manifold, sup- 
pressing its displacements and as a result the manifold 
can be bound to the strain field. Self-localization takes 
place only if the overall free energy of the system with 
the SL manifold is less than that of the free manifold in 
the strain- free medium. 

We start from two examples demonstrating the physics 
of self- localization. 

1. A domain wall in a two dimensional elastic medium. 
The position of the wall relative to axes x and t can be 
represented by the function x — h{t). The interactions 
of the domain wall with the elastic medium can then be 
modeled by a Hamiltonian in which the wall is repre- 
sented as an elastic string of stiffness coefficient m: 

H ^^mj dt{dh/dtf -V J dt[du/dx{x = h)] 

-f^ J dxdt[A{du/dtf + B{du/dxf] (I) 

The deformations of the medium give rise to the last 
term; for simplicity these have been described by a scalar 
displacement field u, with A and B being elastic con- 
stants. The second term provides the coupling between 
the deformation and the wall, in terms of the strain eval- 
uated at the wall position x = h. Since —du/dx is pro- 
portional to the local density change, phenomenologically 
the coupling —V{du/dh) is the lowest order term in the 
expansion of the interaction in powers of the variable part 
of the medium density. The coupling constant V is inde- 
pendent of the elastic properties of the medium, and can 
be estimated in terms of the bulk phase transition tem- 
perature Tc and lattice spacing a a.s \ V\a = Tc. In writing 



the coupling in the form —V{du/dh) it was assumed that 
we are not too close to the bulk critical point so that the 
range of the interaction between the wall and an indi- 
vidual atom is of order a, which is invisibly small in the 
long-wavelength theory we are using and thus does not 
appear; however the scale a will provide a short-distance 
cutoff implicitly present in Eq. (I). 

Placing the wall at the origin h — Q, and minimiz- 
ing Eq. (1), we find du/dx — {V/B)5{x) and u{x) = 
{V/2B)sign{x), which means that at T = the strain 
field is localized at the domain wall (the true range of 
the delta-function is the lattice spacing a). The magni- 
tude of the maximal strain can be estimated as |y|/-Ba = 
\V\a/Ba'^ ~ Tc/Ba^] it is small compared to unity, since 
Tc is much smaller than the binding energy per atom 
Ba^] this justifies the use of linear elasticity theory. 

At a finite temperature T, the wall and the medium 
degrees of freedom each fluctuate strongly - for the un- 
coupled system {V = 0) the mean-square fluctuations 
diverge with the system size L SiS < h? >= TL/m, and 
< u2 [T/VAB] ln(L/a), respectively^. To under- 
stand the system at finite temperature, below we develop 
a variational theory valid to zero order in < > / < 
h'^ > (i.e. ignoring medium fluctuations), and having the 
spirit of Flory-type theories of polymer physics Q . 

Assume that the magnitude of the medium deforma- 
tion is of order uq in a region of size L perpendicular 
to the wall and negligible elsewhere. The correspond- 
ing strain field creates a localizing potential for the wall 
[second term of Eq. (I)]. Both the localizing potential 
and the wall itself fluctuate; however, for a macroscopic 
system the medium fluctuations are much weaker than 
those of the free wall, and so for a first approximation 
we can ignore medium fluctuations altogether, and con- 
sider the domain wall to be fluctuating in a fixed average 
self-consistent potential due to the strain in the medium. 
The localizing potential confines the domain wall within 
a region of typical size L. This diminishes the available 
number of wall configurations, thus decreasing the en- 
tropy and increasing the free energy. The corresponding 
confinement free energy Q is given by T^/mL^, so that 
the total variational free energy (per unit length in the t 
direction) is estimated as 
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where the second term is the free energy gain of localiza- 
tion [corresponding to the second term of Eq. (1), while 
the last term is the increase in the elastic energy due to 
the creation of the strained region near the wall [the last 
term of Eq. (1)]. The variational parameters L and uo 
can be found by minimizing Eq.(2). The minimum over 
Mo is given by uq = V/B, and then Eq.(2) becomes 



F{L,uo = V/B)^ 
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This expression (shown schematically in Figure la) has 
a minimum for finite L and Faun < 0, which implies that 
despite strong thermal fluctuations an elastic strain al- 
ways accompanies the domain wall; the strain is localized 
within an equilibrium localization length given by 



Leg ^ T'^B/mV'^ 



(4) 



We will call the domain wall (plus the accompanying 
strain field) self-localized because the wall is "attached" 

to the strain field it generates - the probability of finding 
the wall a distance x away from the strain peak is signif- 
icant for |a;| < Leg, negligible elsewhere, and goes to zero 
as \x\ — > oo. Note that self-localization docs not imply 
any breaking of translational symmetry; both the wall 
and the strain field can be shifted simultaneously without 
energy cost. Fluctuations of the SL wall will necessarily 
involve the degrees of freedom of the medium; therefore 
the SL wall possesses a larger stiffness than the same free 
wall in an undcformable medium. 

The continuum result (4) is valid whenever the local- 
ization length Leg is much bigger than the lattice spacing 
a. Correspondingly we will distinguish between the cases 
Leg ^ a (the macroscopically SL state) and L^q = a 
{microscopically SL). These meet at the crossover tem- 
perature T* = \V\y/ma/B < T^. For T < T* wc have 
Leg = a. Since the domain wall is an anisotropic object, 
it is also localized in the t direction over a different length 
scale Lt, which can be estimated from L^^ = TLj/m. At 
the crossover temperature T*, we find Lt* = ma^ /T* = 
{Tc/T-*)a ^ a. On length scales smaller than Lt the do- 
main wall is essentially decoupled from the continuum, 
and so the fluctuations of the medium at that scale can 
be estimated as < >= [T/ \J AB\ \n{Lt/a). Our theory 
will be valid whenever < >, evaluated at the smallest 
allowed scale Lt*, is much smaller than a^. This gives 
rise to the inequality 



T* 



: ln{ma/T*) < 1 



(5) 



For most practical purposes {A = B) the condition (5) 
is satisfied exceptionally well, because T* <c Tc ^ Ba^. 
Only when there are special reasons for A to be very small 



(a quasi-one- dimensional medium, for example), do we 
have to start worrying about the role of the medium fluc- 
tuations. For every case consider below the neglect of the 
medium fluctuations can be justified on similar groimds. 

The consequences of manifold self-localization depend 
on the dimensionalities of the manifold and the contin- 
uum as demonstrated by the next example. 

2. A flux line (directed polymer) of stiffness m in a 
three- dimensional crystal. We will assume that the 
strains are localized within a region of size L around the 
vortex. The energy cost of the elastic distortion [the last 
term in Eq.(2)]. is now estimated as BL^{uo/L)^ = Buq 
where B is some combination of the elastic constants of 
the crystal. As a result, the total variational free energy 
per unit length will have the form 



F(L,uo) = 
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The minimum value with respect to uo is 

F{L,uo = V/BL)~' 
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(6) 
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In contrast with Eq. (3), both terms have the same scal- 
ing dependence on L, and the outcome depends on the 

relative size of the coefficients. The two terms are of 
similar size at the critical temperature 



(8) 



for T > Tioc the free energy (7) is minimized for the 
free vortex (L = oo) in the strain- free medium [uq = 0), 
while for T < Tioc the free energy is minimized as L ^ 0. 
Since the strain (which is of order uq/L = V/BL"^) di- 
verges as L ^ 0, this indicates that the linear elasticity 
theory used to estimate the last term of Eq.(6) fails be- 
low the localization temperature (8). There are several 
equivalent ways to cure this problem. 

First, we can add anharmonic terms to the right-hand 
side of Eq. (6) which will guarantee that the equilibrium 
L cannot fall below some microscopic scale for T < Tioc- 
Alternatively we can introduce a short-length cutoff a 
and continue to use Eqs.(6) and (7), but restrict the mini- 
mization to the region L > a. Neither of these procedures 
provides a quantitative picture for T < Tioc but for the 
purposes of illustration we note that using the cutoff (we 
adopt it hereafter) predicts that for T < Tioc the free 
energy density is given by 



F^{T^-T{'oc)/ma' 



(9) 



In contrast to the case of a domain wall in two dimen- 
sions, a flux line in a three-dimensional crystal can exist 
either in the high-temperature free state or in the low- 
temperature microscopically SL phase. The theory de- 
veloped here is valid only for temperatures well below 
the superconductive temperature Tc, since Tc <C Ba^, 
one can deduce from (8) that Tioc ^ Tc. 
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The nature of these results - having the same 
dependence of both the "localizing" and " delocalizing" 
terms of the free energy (7), plus the fact that the expres- 
sion for the phase transition temperature (8) is cutoff- 
independent - shows that three dimensions is the lower 
critical dimensionality of the self-localization problem. 

The critical behavior (9) predicted by our variational 
theory is only an estimate, and a more detailed theory 
is needed to reveal the essential singularities that are ex- 
pected to occur at the critical point in a marginal case. 
The long-ranged nature of bulk elasticity is clearly im- 
portant here: the equilibrium medium displacement away 
from the vortex depends on the variational localization 
length L as uo = V/BL[see (7)], implying that the radial 
component of the displacement field depends inversely 
(u = V/ Br) on the distance r from the vortex. 

3. Now let us consider the general case of a D- 
dimensional manifold of stiffness m placed in a d- 
dimensional continuum. 

For D < 2 the variational free-energy density general- 
izing Eqs.(2) and (6) has the form 
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where the L-dependence of the first term ||] demonstrates 
the general tendency that thermal fluctuations become 
less important with increasing manifold dimensionality 
D. The case D = 2 is the upper critical dimensionality |^] 
for thermal fluctuations and will be treated separately. 
For D > 2 thermal fluctuations are unimportant, and 
the free energy density will be given by just the last two 
terms of (10). Eq.(lO) includes several important special 
cases, such as a vortex line or directed polymer {D = 1), 
and a domain wall {D = d — 1), correspondingly placed 
in a d-dimensional elastic continuum characterized by an 
elastic constant B. 

Eq. (10) is minimized (with respect to uq) for Umin = 
V/ BL"^^^^^ , which immediately implies that only for the 
case of a domain wall {D — d — 1) is the induced strain 
field well-localized; in all other cases it is long-ranged. 
The corresponding minimum value of (10) is 



F(L,Mo 



= 



D/(2-D) 



2BL'i-D 



(11) 



Here we have several cases to consider: 

For D > 2 (or for any D and T = 0) the first term of 
(11) will have to be dropped, and what is left is minimal 
for L ^ 0. We must assume there is a short-range cutoff 
a, and conclude that the manifold will be in a microscop- 
ically SL state with L^q = a. 

For 2D/{2- D) > d - D and D < 2 the right-hand 
side of (11) (shown schematically in Figure la) has a 
minimum at 



[(BT/V2)2-^(T/m)^] 



Dil/[2D-((i-D)(2-_D)] 



(12) 



which generalizes Eq.(4). At the minimum one has 
-Fmin < 0, which implies that the manifold is in a macro- 
scopically SL state as long as L^q S> a. Should the local- 
ization length Lgq decrease (due to some change in the 
system parameters), the system continuously crosses over 
into a microscopically SL state characterized by Leg — i- 

The case 2D/{2-D) = d-D and D <2, similar to the 
previously analyzed problem of a flux line, is marginal. 
Here the manifold can exist either in a high-temperature 
free state or in a low- temperature microscopically SL 
phase separated by a continuous phase transition. 

For 2D/{2~ D) < d~ D and D <2 the right-hand 
side of (11) (shown schematically in Figure lb) has a 
maximum at a scale L^ax [actually given by (12)], ap- 
proaches zero from above as L — > oo, vanishes at a scale 
Lo < imax[up to a scalc factor Lq is still given by the 
right-hand side of (12)], and goes to minus infinity as 
L ^ 0. Here the outcome depends on the relationship 
between the cutoff a, and the scales Lq and Lmax- 

For a > Lmax (the high-temperature limit) the free- 
energy density (with the cutoff) has a single minimum at 
L ~ oo: the manifold is in a free state. 

For Lq < a < Lmax (intermediate temperatures) the 
free energy density has a global minimum at L — oo, and 
a local minimum ai L ^ a. The ground state is still a 
free manifold, but now a microscopically SL state can be 
metastable because of the barrier separating it from the 
free state. 

For a < Lq (a low-temperature limit) the free energy 
density has a global minimum at L = a, and a local min- 
imum at L = oo. The ground state is a microscopically 
SL manifold but a free manifold can be metastable due 
to the free energy barrier separating it from the SL state. 

The condition a = Lq determines a first-order phase 
transition point between the microscopically SL phase 
and a free manifold. 

4- For the marginal case D = 2 we will consider only 
the practically most important case of a domain wall in 
a three- dimensional crystal. The free energy per unit 
area will now have the form 

F(L, uo = V/B) = (T/a2)e-'"^'/^ - V^/BL (13) 

where the first term is the free energy of confinement Q] 
which for D = 2 replaces the first terms of Eqs. (10-11). 

To simplify the analysis it is convenient to in- 
troduce the dimensionless coupling constant 7 — 
{V'^a^/B){m/T'^Y/'^. For D = 2 the parameter V is es- 
timated as \V\a'^ = Tc = ma^, and the corresponding 
couphng constant is 7 = {Tc/Ba^){Tc/T)^/^ , thus imply- 
ing that for T < Tc the parameter 7 can vary from a very 
small value of order Tc/ Ba^ <C 1 (for a temperature close 
to Tc) to infinity (the zero-temperature limit). 

In the low-temperature limit (7 ^ 1) the func- 
tion (13) is monotonic increasing which means that the 
ground state is a microscopically SL domain wall. In 
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the high temperature hmit (7 <C 1) the right-hand 
side of (13) (shown schematically in Figure Ic) has 
a maximum at some Lmax, and a minimum at some 



> in 



The transition between the two behav- 



iors occurs at 7 = 1 which corresponds to the temper- 
ature Ti ^ Tc{Tc/Ba^)^/'^ < Tc and Li = L^ax = 
L,„i„ ^ (ri/m)i/2. Since the scale Li ^ {Ti/my/^ ^ 
a(Ti/Tc)^/^ is much smaller than the cutoff a we con- 
clude that the domain wall will be in a microscopi- 
cally SL state L^q = a both for T < Ti and in some 
temperature range above Ti until the cutoff a coin- 
cides for the first time with iinin(^min increases with 
the temperature). This will happen at some 7 <C 1, 
and in this limit Lmin can be estimated as Lmin — 
(T/to)^/^ ln^/^(l/7). The condition Lmin = a determines 
a crossover temperature T* as a solution to the equa- 
tion ln"'^/^(l/7) ~ a(m/T)^/^. Remembering that 7 — 
{V^ayB){m/TY^^ = {TJ Ba^){TjTf/^ , we can solve 
the equation with logarithmic accuracy, \u{Ba?' /Tc) S> 1, 
to find T* ^ Tc/lTi{Ba^/Tc). 




FIG. 1. Sketch of the free energy F of the manifold as a 
function of the variational localization length L for the var- 
ious regimes of manifold and medium dimensionalities: (a) 
2D/{2-D) > d-D, D <2; (b) 2D/{2-D) < d-D, D < 2; 
(c) a domain wall in a three-dimensional crystal (D=2,d=3) 
in the high-temperature limit. 

The physical significance of the crossover temperature 
T* is that for T < T* the free energy density has a single 
minimum at L^q = a, and the domain wall is in a micro- 
scopically SL state. At the same time for T > T* the 
free energy density has a minimum at 



^ (T/m) \n{B^TymV^a'^) 



(14) 



The ground state is still self-localized but its localization 
length (14) is bigger than a and grows with tempera- 
ture. The bulk phase transition temperature Tc sets a 
natural upper bound to Leg, which can be estimated as 
Leq < a\n^^^{Ba^ /Tc). In reality the dependence (14) is 
only reliable for T* < T < Tc excluding the temperature 
range of bulk critical fluctuations. 

Because thermal fluctuations are only marginally rel- 
evant for a domain wall in a three-dimensional crystal, 



experimental verification of the law (14) (for example, 
by measuring the strain profile induced by the wall) may 
not be easy unless \n{Ba^ /Tc) » 1, and the bulk critical 
region is narrow. 

It is known experimentally and well-understood the- 
oretically [H how topological defects and other inhomo- 
geneities deform the embedding medium at zero temper- 
ature. In our language this corresponds to the microscop- 
ically SL manifold. As we have shown, thermal fluctu- 
ations can change this picture drastically leading to the 
possibility of macroscopically SL states or even strain- 
free states. The fact that the manifold is in a SL state 
should be especially important when the accompanied 
strain field is long-ranged (this excludes domain walls), 
and when there are many manifolds present. We specu- 
late that the SL transition predicted for a single flux line 
will manifest itself in a structural transition of a dilute 
system of vortices. The present theory implies that there 
is a similar phenomenon for any D = 1 manifold in a 
deformable medium in three dimensions. 

The problems analyzed above involve directed mani- 
folds placed in a deformable medium. For the special 
case of line manifolds {D = 1) they are connected to 
quantum-mechanical phenomena at T through the 
Feynman formulation of quantum mechanics The 
relationship goes through the correspondence tempera- 
ture — > Planck's constant, line configuration — s- particle 
world line, position along the line — > imaginary time, and 
line stiffness particle mass. Then line self-localization 
phenomena are formally related to the polaron effect in 
quantum physics - an electron placed in an elastic con- 
tinuum might get self-trapped by its own strain field 
More precisely, the Hamiltonian Eq. (1) can be viewed 
as an imaginary-time Action for the polaron problem in 
one dimension, while the expression for the localization 
length (4) can be related to the size of a large polaron in 
the adiabatic limit 0. This implies that the SL mani- 
fold has a well-localized probability distribution function 
which depends only on transverse coordinates. This is 
qualitatively different from the free manifold, which has a 
Gaussian probability distribution function that depends 
on all coordinates. 
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